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An epoch of Higgs relaxation may occur in the early universe during or immediately following 
postinfiationary reheating. It has recently been pointed out that leptogenesis may occur in minimal 
extensions of the Standard Model during this epoch [1]. We analyse Higgs relaxation taking into 
account the effects of perturbative and non-perturbative decays of the Higgs condensate, and we 
present a detailed derivation of the relevant kinetic equations and of the relevant particle interaction 
cross sections. We identify the parameter space in which a sufficiently large asymmetry is generated. 


I. INTRODUCTION 

During the inflationary era, the Higgs field may de¬ 
velop a stochastic distribution of vacuum expectation 
values (VEVs) due to the flatness of its potential [2-9], 
or it may be trapped in a quasi-stable minimum. In 
both cases, the Higgs field relaxes to its vacuum state 
via a coherent motion, during which time the Sakharov 
conditions [10], necessary for baryogenesis, are satisfied 
by the time-dependent Higgs condensate and the lepton- 
number-violating Majorana masses in the neutrino sec¬ 
tor. At large VEVs, the Higgs field may be sensitive to 
physics beyond the Standard Model, which can generate 
an effective chemical potential which increases the energy 
of antileptons in comparison to leptons. In Ref. [1], we 
used a Oq operator familiar from spontaneous baryoge¬ 
nesis models (e.g., [11]) to produce a baryon asymmetry 
matching cosmological observations. 

In this work, we build on our previous analysis. In 
particular, we replace the estimate of the Higgs-neutrino 
cross section with a tree-level calculation which includes 
resonant effects. Additionally, we include the effects of 
Higgs condensate decay, with both perturbative and non- 
perturbative contributions. We also present a detailed 
derivation of the relevant Boltzmann equation. Addition¬ 
ally, the effective Oq operator can be generated through 
fermionic loops; therefore, its scale can be set either by 
some heavy mass scale or by the temperature of the 
plasma. We consider additional combinations of these 
scales along with mechanisms to generate the large Higgs 
VEV during inflation, and in particular, we also present 
an analysis of the relevant parameter space. 

We note that as in Ref. [1], we consider an asymme¬ 
try produced via the scattering of neutrinos and Higgs 
bosons in the plasma produced by the decays of the in- 
flaton, and therefore this scenario requires a relatively 
fast reheating. This is in contrast to Ref. [12], which 
similarly considered the same Oq operator but produced 
the matter asymmetry via the decay of the Higgs con¬ 
densate. 

While we focus here on the relaxation of the Higgs 
field, it has also been observed that the axion field can 


undergo a similar post-inflationary relaxation [13, 14], 
and our analysis can easily be extended to this scenario. 

The structure of this paper is as follows. In the next 
section, we consider two specific mechanisms by which 
the Higgs field can acquire a large vacuum expectation 
value during inflation. In the scenarios considered here, 
the subsequent evolution of the Higgs VEV produces an 
effective chemical potential, which influences the inter¬ 
actions of leptons in the thermal plasma produced via 
reheating. The presence of the plasma, however, also 
influences the evolution of the Higgs VEV through fi¬ 
nite temperature corrections to the effective potential. 
Therefore, we discuss reheating in section HI before we 
consider the evolution of the Higgs condensate in section 
IV. Next, we introduce a higher-dimensional operator, 
involving only Standard Model fields, which represents 
new physics at some high energy scale. In section V, we 
demonstrate that, while the Higgs VEV is in motion, this 
operator induces an effective chemical potential which 
distinguishes leptons from antileptons. We derive the 
resulting Boltzmann equation for lepton number in sec¬ 
tion VH. In section VHI, we present a numerical analysis 
covering a variety of initial conditions and scales for new 
physics, and we identify the allowed parameter space for 
a successful leptogenesis. 

II. INITIAL CONDITIONS FOR THE HIGGS 
VEV 

We begin by motivating our project with the obser¬ 
vation that the Higgs field can acquire a large vacuum 
expectation value (VEV) for a variety of reasons during 
inflation; therefore, an epoch of post-inflationary Higgs 
relaxation is a general feature of many cosmological sce¬ 
narios. In this work we are interested in generating an 
excess of leptons over antileptons during this epoch. 

When we generate the lepton asymmetry during this 
epoch of Higgs relaxation, we will find that the result¬ 
ing asymmetry depends on the initial value of the VEV, 
denoted yj (0^) = 4>q. During inflation, quantum fluc¬ 
tuations of the Higgs field were ongoing, and therefore 
different patches of the Universe had slightly different 
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VEVs at the end of inflation. Regions that begin with 
slightly different (j)Q values consequently develop different 
baryon asymmetries. This produces unacceptably large 
baryonic isocurvature perturbations [15], which are con¬ 
strained by CMB observations [16]. Therefore, in order 
to suppress isocurvature perturbations in the late uni¬ 
verse, it is necessary to have a small variation in these 
values. 

In this section, we discuss two ways of generating the 
requisite large VEVs while suppressing the variation be¬ 
tween different spacetime regions: through quantum fluc¬ 
tuations, which are suppressed by a Higgs-inflaton cou¬ 
pling until the end of inflation, and by trapping the Higgs 
field in a false vacuum. 

The Standard Model Higgs boson has a tree-level po¬ 
tential 


V($) (1) 

where $ is an SU(2) doublet. The classical field may be 
written as 
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where (^(cc) is a real scalar field. The parameters m and A, 
although constant at tree-level, are modified by both loop 
and finite temperature corrections. For the experimen¬ 
tally preferred top quark mass and Higgs boson mass, 
loop corrections result in a negative running coupling 
A at sufficiently large VEVs, with the result that the 
(j) = vew = 246 GeV minimum is metastable at zero tem¬ 
perature [17]. We note, however, that a stable vacuum is 
possible within current experimental uncertainties [17]. 

The running of the quartic coupling produces a shallow 
potential, which has the consequence that that a large 
VEV develops during inflation due to quantum fluctua¬ 
tions, at least in the regime in which the Standard Model 
vacuum is stable [8] . We consider this sort of scenario in 
subsection IC-2 below. Alternatively, the metastability 
of the electroweak vacuum is frequently possible within 
the inflationary paradigm [18], and the Higgs potential 
may be sensitive to higher-dimensional operators which 
lift the second minimum. We consider this scenario in 
the subsequent subsection. 


A. IC-1: Metastable Vacuum at Large VEVs 

At the large VEVs, the Higgs potential may be sensi¬ 
tive to the effects of higher-dimensional operators, which 
can lift the second minimum and consequently stabilize 
the electroweak vacuum. The Higgs VEV may take an 
initial large value during inflation, similar to the initial 
VEV of the inflaton field itself in chaotic inflation mod¬ 
els. During inflation, such a VEV evolves towards the 
false vacuum from above, and then remains trapped in 
this false vacuum until destabilized by thermal correc¬ 
tions in reheating. Subsequently, the field rolls to the 


global minimum at (/) = 0, until electroweak symmetry is 
broken at a significantly later time. 

In order to lift the second minimum, we consider terms 
of the form 


“ A6 ■ 

This non-renormalizable operator may be viewed as an 
effective operator arising from integrating out heavy 
states in loops. 

During inflation, thermal corrections in the su¬ 
percooled universe are insufficient to destabilize the 
metastable vacuum. We also ensure that the quantum 
fluctuations (discussed in detail in the next subsection) 
do not destabilize the vacuum by requiring that the po¬ 
tential barrier height AV ^ Hj. In order to suppress the 
above-mentioned isocurvature perturbations, we will en¬ 
sure that fluctuations about the false minimum are able 
to relax back to the minimum, for which it is sufficient 
to ensure iries y/cPVjd^ > Hj in the region probed 
by quantum fluctuations. 

As a specific example, we consider the Higgs potential 
with one loop corrections [17] with the experimentally 
preferred values mh = 126 GeV and mt = 173.07 GeV. 
Taking Auft = 6.52 x 10^® GeV gives a metastable mini¬ 
mum near (p = 10^^ GeV, with a potential barrier height 
of AV « 10®^ GeV^. We will consider Hi ~ 10^^ GeV; 
in addition to being insufficient to probe the region be¬ 
yond the barrier, this is less than the effective mass 
rueff 10^^ GeV in the region probed by quantum fluctu¬ 
ations. Provided that the maximum reheat temperature 
is greater than ~ 5 x 10^^ GeV, thermal corrections dur¬ 
ing reheating are sufficient to destabilize this vacuum. 


B. IC-2: Quantum Fluctuations 

The running coupling constant A results in a shallow 
potential, and during inflation, scalar fields with slowly 
rising potentials generically develop large VEVs. Quali¬ 
tatively, the scalar field in a de Sitter space can develop 
a large VEV via quantum effects, such as Hawking-Moss 
instantons [2, 3] or stochastic growth [4-6]. The field 
then relaxes to its equilibrium value via a classical mo¬ 
tion, which requires a time 

T0 ~ m~g ~ . (4) 

If the universe expands sufficiently quickly during in¬ 
flation, then relaxation is too slow and quantum jumps 
occur frequently enough to maintain a large VEV. Specif- 
ic ally, l arge VEVs occur if the Hubble parameter Hj = 
y^87r/3Aj/Mpi ^ For field values p that satisfy 

this relation, Hubble friction is sufficient to prevent the 
system from relaxing to its equilibrium value p = 0. Av¬ 
eraged over superhorizon scales, the mean Higgs VEV is 
such that V^pi) ~ Hj [2, 3, 8], provided that this VEV 





3 


does not probe the second vacuum in the case that the 
electroweak vacuum is quasistable. 

Although the average vacuum expectation value is (pi, 
there is variation between the VEVs of different horizon¬ 
sized patches. Consequently, different patches of the ob¬ 
servable universe began with different cpo values, and as 
discussed above, this generically results in unacceptably 
large isocurvature perturbations. However, also as men¬ 
tioned, the Higgs potential is sensitive to the effects of 
higher-dimensional operators at large VEVs; here we use 
such operators to limit the growth of the Higgs VEV to 
the last several e-folds of inflation. This has the result 
that the isocurvature perturbations are limited to smaller 
angular resolution scales than have been experimentally 
probed. Specifically, we introduce one or more couplings 
between the Higgs and inflaton field of the form 

j^^m+n-2 > 

which increases the effective mass of the Higgs field dur¬ 
ing the early stages of inflation, when (/) is large (super- 
planckian, in the case of chaotic inflation). As explained 
above, when ~ ~ Hj the expansion of the 

universe is not sufficiently rapid to trap the field at large 
VEVs. At the end of slow-roll inflation, (/) decreases; 
consequently this term becomes negligible and the Higgs 
acquires a large vacuum expectation value. 

If the Higgs VEV grows during the last Mast e-folds of 
inflation, it reaches the average value 

(po = min[M, \/MastM/27r]. (6) 


Provided Mast « 5 —8, the baryonic isocurvature pertur¬ 
bations develop only on the smallest angular scales which 
are not yet constrained. 

We emphasize that operators of the form (5) may be 
viewed as effective operators arising from integrating out 
heavy states in loops. We note that the change in (/) dur¬ 
ing the slow-roll phase of inflation is model-dependent, 
and consequently the allowed range of parameters c, m, 
and n differs from model to model. This range may be 
quite narrow, and so this scenario may require some fine- 
tuning. 

As a concrete example, we consider only the term. 
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which induces an effective mass meff((/)) = (/)" 
for the Higgs field. We define Ii as the VEV of the in¬ 
flaton field value at the end of slow roll inflation, and 
I 2 as the VEV of the inflaton field 8 e-folds before the 
end of slow roll inflation. To ensure that the Higgs VEV 
grows only during the last e-folds, we must choose pa¬ 
rameters such that rriesih) ~ Hj. We illustrate this ap¬ 
proach with quartic inflation (although this is disfavored 
observationally; see Ref. [16]). With the inflaton poten¬ 
tial Vi = A//"*, slow roll inflation ends when the inflaton 


as a vacuum expectation value of Ii = Mpify/^. The 
number of e-folds during the time in which the inflaton 
evolves from (J) to Ii is 

which gives 



(Although this field value is superplanckian, this is a fea¬ 
ture of quartic inflation which does not necessarily apply 
to other inflationary models.) The Hubble parameter at 
this field value is given by 

The quartic coupling A/ must be < 10“^^ in order to 
avoid large CMB temperature anisotropies, which gives 
M ~ lO^Mpi for both n = 2 and n = 4. In this way, a 
coupling between the Higgs field and the inflaton field can 
prevent the Higgs VEV from growing until the last sev¬ 
eral V-folds of inflation, suppressing the scale of isocur¬ 
vature perturbations. As this example illustrates, the 
constraints on the Higgs-inflaton coupling depends on 
the shape of the inflaton potential. Although we have 
demonstrated an explicit calculation using a quartic in¬ 
flationary potential, a similar calculation can be done 
with other potentials. 

Our analysis of the final asymmetry will depend only 
on the VEV of the Higgs field at the end of inflation, 
which as noted above is (po = min]^/, -v/MastM/STr] pro¬ 
vided that the parameters in whatever inflationary model 
is used are chosen such that the Higgs VEV does not be¬ 
gin to grow until the last Mast-folds of inflation. There¬ 
fore, we do not specify a specific inflationary model in 
our analysis, and we take the inflationary scale A„ and 
the inflaton decay rate T/ to be free parameters. 


III. REHEATING 

Now that we have established that the Higgs field can 
develop a large VEV during inflation, we are interested in 
its subsequent evolution to its equilibrium value. Relax¬ 
ation begins when the Hubble parameter is comparable 
to the effective mass of the Higgs field, mef[{(p) « H{I), 
which is within the reheating epoch. Therefore, this re¬ 
laxation is sensitive to finite temperature effects due to 
the plasma, and so we now proceed to discuss reheating. 

In both scenarios, whether IC-I or IC-2, we ensure that 
the energy density is never dominated by the Higgs field. 
Inflaton oscillations dominate until the transition to the 
radiation dominated era, which occurs when the infla¬ 
ton decay width is comparable to the Hubble parameter, 
E/ ~ Hrh, which typically occurs after the Higgs field 
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has lost a significant portion of its energy. Consequently, 
the reheat temperature Trh ~ vTTA^r, is generally only 
weakly constrained [19]. 

For simplicity, we assume coherent oscillations begin 
instantly at the end of the inflationary epoch, and as a 
simple model, we assume that the inflaton decays entirely 
to radiation at a constant rate. 


until it reaches the reheat temperature Tr ^ \/TiMp\. 
Subsequently, radiation dominates the energy density 
and the temperature evolves as 



Pr + 4:H{t)pr = Tipi, 


IV. EVOLUTION OF THE HIGGS VEV 


where 


We now turn our attention to the relaxation of the 
Higgs VEV, which evolves as [20] 


PI 




( 12 ) 


is the energy density of the inflaton field. The evolution 
of the Hubble parameter is given by 


^ I='j +"'I 

This is a complete system of equations that may be solved 
independently of the evolution of the Higgs condensate. 
Throughout this work, we take t = 0 to be the beginning 
of the coherent oscillation of the inflaton field; during 
the coherent oscillation epoch, the universe evolves as if 
it were matter dominated, until the radiation from re¬ 
heating dominates. 

During reheating, the effective temperature of the 
plasma is defined using the radiation density as 

9^11 rj-i4 /I 

= (14) 

For t ^ ti = (2/3) -^/S/SttMpi /A^, the temperature 
evolves as 


Vg*7r3 t ) 



+ + = (17) 

where V^{(j),T) is the Higgs effective potential, includ¬ 
ing modifications from the decays of the condensate [8]. 
F// describes the effect of the perturbative decay of the 
condensate. In the first subsection below, we discuss the 
one-loop corrected potential, including one-loop correc¬ 
tions to the RG equations. Subsequently, we consider the 
non-perturbative decay of the Higgs condensate, followed 
by perturbative decay. Finally, we present a numerical 
analysis of the evolution of the Higgs condensate, before 
we proceed to the next section which introduces the rel¬ 
evant higher dimensional operator we use to produce the 
nonzero lepton asymmetry. 


A. Effective Potential 

The Standard Model Higgs potential computed to a 
fixed order in perturbation theory is generally gauge- 
dependent, although the value of the potential at the 
extrema are not (see, for example, [21, 22]). One can 
ensure gauge-invariant results by removing the gauge- 
dependence of the potential using Nielsen identities [23- 
25]. Here we use the Landau gauge, which has good nu¬ 
merical agreement with the corrected potential [21, 26]. 
In our analysis, we have used the one-loop corrected po¬ 
tential [27], with running couplings (including one-loop 
corrections to the renormalization group (RG) equations, 
as given in [17]). The one-loop potential is 


T 7-1—loop 2 i2 


4 _|__ 

'4"^ (47r)2 


(rnnicji)'^ 


4 


In 


3 

’ 2 
V 


3tog(</>)'‘ ^ wg(^)^ 




(18) 


where p, is the renormalization scale and the tree-level 
masses for the Higgs boson, Goldstone mode, W bosons. 


Z boson, and top quark are 




4 

2 




mt = 


mH = mi+3\4>, mQ=mi + \. 


72’ 

(19) 
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We have also included the finite temperature correc¬ 
tions [28, 29], 


where the kinetic terms of the gauge fields can be ex¬ 
panded as 




T2 

27r^ - 




(^) 




+\2mpF (^)j 



where 


LXJ - 

JB{y) = -^K2{ny), ( 21 ) 

n—l 

My) = E ^2(ny), (22) 


Ca, kin = - ^ (V^w+ - V,W+) (V'^W - V"W"-) 

1 2 

— - {Vfj,Ziy — ViyZf^) + 0 (g)(non-Abelian terms). 

(24) 

Since the non-Abelian contributions are small at the be¬ 
ginning of the resonant production of W and Z bosons, 
we ignore those terms [ 8 , 30]. We also work specifically in 
flat FLRW spacetime, 5 ^ 1 / = {T)r]f^^, with conformal 

time T = f a~^dt. The resonant production of the weak 
gauge fields, 01 Z^, in momentum space is then 

described by 


and we have ignored the contributions from Higgs bosons 
and Goldstone mode, which only dominate when (j) ^ 
uew- We emphasize that we do not use the high tem¬ 
perature expansion, as during reheating the condition 
T{t) ^ is not satisfied at all times. The renor¬ 
malization scale fi is taken to be \/(p -I- T^. 

We note that two-loop corrections may be significant 
at the boundary of the metastability region [17]; however, 
a self-consistent analysis at two-loop order would include 
finite temperature effects in the RG equations, which is 
beyond the scope of this work. 

After the Higgs VEV passes through zero, it generally 
oscillates around its minimum at ^ = 0 , which remains 
a minimum for T 3> uew- During this oscillation, the 
Higgs condensate can then decay perturbatively and non- 
perturbatively into Standard Model particles. The non- 
perturbative decay happens much faster than the per¬ 
turbative decay and is the dominant channel, as pointed 
out by [ 8 ]. We now proceed to discuss the effect of these 
decays. 


B. Non-Perturbative Decay 


First, we consider non-perturbative decay. The oscil¬ 
lation of the Higgs field provides a time-dependent mass 
term for all the coupled particles, which can cause reso¬ 
nant production of the particles. The produced particles 
then induce an effective mass term to the Higgs conden¬ 
sate as a backreaction; this attenuates the oscillation of 
the Higgs field until the resonant production is off [ 8 , 30]. 

The non-perturbative decay channel of Higgs is domi¬ 
nated hy h —>■ WW, ZZ. The Lagrangian containing the 
Standard Model weak gauge fields and the Higgs sector 
is 






+ 0 + g''^) z^z, 


+ 


kin 5 


(23) 


and 


-ikA'^(^k,T^ 


(</>) ’ 

(25) 

{4>) At^ = 0 , 

(26) 


(</*) + 


2k^ 


dr In [am a) = 0 , 


(27) 


where = \/k‘^ + {(f)) and prime denotes differen¬ 

tiation with respect to conformal time dr. At {k, and 

At (^k, are the transverse and longitudinal components 

of the spatial component A ^k, , respectively. The mass 

term m\ (4>) is given in Eq. (19) and can include the ther¬ 
mal correction by replacing cjp ^ (jp CaT^ where we 
use Cw = 2/3, and Cz < 1 is determined by diagonaliz¬ 
ing the mass matrix [31]. Due to extra friction term in 
Eq. (27) for the longitudinal component At, we expect 
the resonance production of this mode to be suppressed. 
Aq, which depends only on through Eq. (25), should 
also be suppressed [30]. Hence, we focus on the trans¬ 
verse mode At only. 

Resonant production of particles can be understood as 
the amplification of vacuum fluctuations. The number of 
particles in each mode produced from the vacuum is 




■ 


At,{k,T) 



(28) 


The initial conditions are taken to be the WKB approx¬ 
imation of the vacuum solution. 


AT{k,0) — A'T{k,Q) — ( 29 ) 

which satisfy nfe(O) = 0 and the Wronskian condition 
AA'* - A*A' = i. 

Fig. 1 shows the amplification of the W field, which 
increases each time the Higgs VEV passes through zero. 
















6 



4 r 


FIG. 1. Real and imaginary parts (blue and purple lines) of 
Wt{t) for fc = 0 for IC-1, with the parameters A/ = 10^® GeV 
and F/ = 10® GeV. The vertical lines designate the first 
time the Higgs VEV crosses zero, and the time of maximum 
reheating, from left to right. 



FIG. 2. log [nfc(T)] for fc = 0, with the same parameters as 
Fig. 1. Note Uk (r) stops increasing at r ~ 3300/<)>o because 
the amplitude of (j) has decreased to the order of T. The 
effective mass of W was then dominated by T instead of 4>. 


The number density Uk is shown in Fig. 2. It has a se¬ 
quence of flat steps, which are separated by peaks. Those 
peaks occur when (/> = 0 ; due to the rapidly chang¬ 
ing mass, the number of particles is not well defined at 
these points. Particle number is well defined only when 
(j) reaches a local maxima or minima. We approximate 
the particle number of quanta within each oscillation 
of (p by its value when (/) = 0 , which is supported by 
the flatness of the steps in Fig. 2. The resonant pro¬ 
duction begins once the Higgs VEV starts to oscillate at 
T ~ 800/())o. The decrease of nfc at t ~ 2500/())o indi¬ 
cates the system has a stochastic resonance, which is a 
distinctive feature of parametric resonance in an expand¬ 
ing universe [32]. The resonant production then ceases 
at T ^ 3300/(()o because the amplitude of (j) has decreased 
to the order of T. 

If we approximate the oscillation of the Higgs VEV by 


4>{t) = 4>m cos(a; 0 r), we can write Eq. (26) as a Mathieu 
equation of the form 

J2 A 

^ ^ + {w? + cos^ z) At = 0 (30) 

where z = m? k. -|- a^m\ (</) = 0,T)] and 
6 ^ « a?m\ = 0) /w^. The Mathieu equation has 

an instability only when b ^ m^. Thus, resonant produc¬ 
tion is suppressed for 

k > fcmax « \/auj4,mA (<^m, 0) - a^m\ (0, T). (31) 

The produced W and Z fields induce an effective mass 
for the Higgs field as a backreaction, 

= (32) 

= —^ ( 9 ^ + 9 '^) {ZfiZ^) (33) 

where the expectation value oi A = W, Z can be approx¬ 
imated as [ 8 , 32] 

, A ^ s -2 / -1 k^dk 

9 " {A,A,) = I —( 34 ) 

In general, the integral in Eq. (34) will need to be reg¬ 
ularized. However, in our case there is no significant 
contribution from Uk values with k ^ kmax-, and so the 
integral is finite. The upper limit can be approximated 
by kmax- One can then include the non-perturbative de¬ 
cay of the Higgs by adding the induced mass terms Eqs. 
(32) and (33) into the Higgs potential in Eq. (17). 

Fig. 3 shows an example of the Higgs evolution with 
the non-perturbative decay for the IC-1 scenario. The 
increasing effective masses from W and Z affect the os¬ 
cillation of Higgs when ^ T; these decrease the 

amplitude of the Higgs oscillation. When the Higgs VEV 
decreases to (/> < T, the resonant production of W and 
Z end, because the non-perturbative decay channel is 
blocked by the large W and Z thermal masses. In this 
case, one has only to consider perturbative decay chan¬ 
nels, discussed in subsection IV C. 

Note the generated W and Z bosons can decay per- 
turbatively into fermions. This decay could in principle 
obstruct the resonant production of W and Z in the usual 
Standard Model case [33, 34]. However, in the parame¬ 
ter space that we are interested in, the average decay 
times of W and Z bosons (riv^z)~^ are longer than the 
semiperiod of the Higgs oscillation. Thus, we have ig¬ 
nored the decay of W and Z in our analysis. 

The analysis in this subsection can be improved by us¬ 
ing a lattice gauge theory [34, 35]. However, as Fig. 3 
demonstrates, the non-perturbative decay of the Higgs 
condensate is relevant only after several oscillations, 
whereas in our scenario the asymmetry will be gener¬ 
ated primarily during the initial oscillation of the Higgs 
VEV. 
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FIG. 3. Non-perturbative decay of the Higgs condensate for 
IC-1, with the parameters Aj = 10^® GeV and Tj = 10® GeV. 
The purple (blue) line corresponds to evolution of the Higgs 
VEV with (without) non-perturbative decay. The brown line 
corresponds to the temperature of the plasma. The vertical 
dashed line indicates the time of maximum reheating. 


C. Perturbative Decay - Thermalization 


which is a convenient form for numerical purposes. We 
will specify the necessary coefficient gr below. In these 
equations, we have made the approximation that the left- 
and right-handed fermions have the same thermal mass 
m (T) = gxT ; generically, this is not true because they 
are in different representations of the Standard Model 
gauge group. However, this difference, which is much 
smaller than the difference between the particle and hole 
contributions, is negligible [31]. 

The dominant fermionic contribution to the thermal¬ 
ization of the Higgs condensate is from the top quark, 
due to the large top-Higgs Yukawa coupling. The ther¬ 
mal mass of the left-handed top quark is [31] 


1 , 3M^ + i(M|-M^) + M2+M2 

(39) 

where Mi are the physical masses at T = 0, and the 
strong coupling gg = 1.220. 

The presence of particles and holes in the fermionic 
plasma provides two thermalization processes for the 
Higgs condensate. A Higgs boson can decay into a pair 
of particles or a pair of holes respectively if 


The perturbative decay of the Higgs condensate is de¬ 
scribed by the friction term Th4 ’ i^i tbe equation of mo¬ 
tion (17). The decay width can be computed through the 
imaginary part of the self-energy operator 


ImH 

^ H = -, 

rriee 


(35) 


rrieff = 2uji (fci); i= p,h (40) 


is satisfied. The contribution of each process to the decay 
width is 


Imlldec 





where mes = Re-^/ {(j), T) /d4>^ is the effective mass 
of the Higgs boson. In a finite-temperature thermal back¬ 
ground, r H corresponds to the thermalization rate of the 
Higgs condensate. 

Here we consider the fermionic decay channels, mo¬ 
tivated by the large top Yukawa coupling. (The domi¬ 
nant bosonic channels, WW and ZZ, are included in the 
non-perturbative calculation, which dominates their per¬ 
turbative contribution.) In the thermal bath of fermions, 
there are additional excitations which are the removals of 
antiparticles from the Fermi sea (holes). The dispersion 
relations for particles and holes are [31, 36, 37] 


9t 

9t 

/l_ 

V k ) 

In 

ujrp k 

k 

2k 

(jjp — k 

9t 

9t 


In 

uih + k 

k 

2k 

V k) 

LOh- k 


= 0, (36) 
= 0, (37) 


where w (k) = LojT, k = k/T, and the subscripts p and 
h refer to particles and holes respectively. Eq. (37) can 
also be expressed as 


ujh = k coth 



tdh + kj 


(38) 


where yt is the top Yukawa coupling, and 


'^h,p 


1 

exp itOh,p) + 1 


(42) 


are the fermion distribution functions. Although this de¬ 
cay channel is blocked when < 2 min [wh (fc)], a Higgs 
boson can also be absorbed by a hole to produce a par¬ 
ticle. The contribution of this absorption channel to the 
width is 


ImH 


abs 


yt 


Y,kl(c:l-kl) [tbl-kj) {UH-Up) 

(43) 


T2 2TTg^ 

where the index i sums over the solutions of 


Weff + W/, {ki) = cOp (ki). (44) 

The total thermalization rate is then the sum of two chan¬ 
nels ImH = ImHabs + InrHdec ■ 

For IC-1, Fig. 4 shows the thermalization rate of the 
Higgs condensate through the top quark compared with 
the Hubble parameter. We see the thermalization rate 
is comparable to the Hubble parameter only after the 
maximum reheating has been reached. Therefore, the 
evolution of the Higgs VEV is affected only at the end of 
reheating later time of reheating as shown in Fig. 4. 
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FIG. 4. Higgs thermalization rate through the top quark com¬ 
pared with the Hubble parameter for IC-1, with the param¬ 
eters A/ = 10^® GeV and F/ = 10® GeV. The vertical lines 
designate the first time the Higgs VEV crosses zero, and the 
time of maximum reheating, from left to right. 



FIG. 5. The evolution of the Higgs VEV for IG-1 (blue line) 
and temperature (purple line) as a function of time, with the 
parameters Aj = 10^® GeV and F/ = 10® GeV. This plot 
includes both the effect of non-perturbative decay and ther¬ 
malization. The vertical line designate the time of maximum 
reheating. 

We have repeated the above analysis with the bot¬ 
tom quark in place of the top quark and verified numeri¬ 
cally that its contribution is negligible; we also note that 
plasma effects can delay thermalization [38]. We also re¬ 
mark that particularly for IC-2, the thermalization rate 
is frequently much smaller than the Hubble parameter. 
Since thermalization occurs on such long time scales, it 
has no effect on our analysis. 


D. Numerical Results 

Figures 5 and 6 illustrate the evolution of the Higgs 
VEV (and temperature) as functions of time. For IC- 
2, the relevant inflaton parameters are A/ = 10^^ GeV 
and F/ = 10® GeV, and we have assumed A^iast = 8 
to determine (j)o, which does not probe the quasistable 
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FIG. 6. The evolution of the Higgs VEV for IG-2 (blue line) 
and temperature (purple line) as a function of time, with the 
parameters Aj = 10^^ GeV, F/ = 10® GeV, and Ni^st = 8. 
This plot includes both the effect of non-perturbative decay 
and thermalization, although the effect of condensate decay 
is not appreciable in this case. 


vacuum. For IG-1, we have used the operator (with the 
numerical parameters) discussed in section H A to lift the 
second minimum, along with the inflationary parameters 
AI = 10^® GeV and F/ = 10® GeV. For both plots, we 
use 126 GeV and 173.07 GeV for the masses of the Higgs 
boson and top quark respectively. 

We briefly remark on the qualitative features of these 
plots. Although in IC-2, the Higgs VEV is constrained 
to grow only in the last 8 Wfolds of inflation, it may still 
reach a large value if Hj is large; for these inflationary 
parameters, Hj = 2 x 10^® GeV. For other choices of A/, 
the initial VEV (j)o for the IC-2 scenarios is significantly 
smaller. Conversely, in our IC-I scenario, the initial VEV 
(j)o « 10^® GeV is set by the parameters chosen in section 
HA. 

In scenario IC-I, the Higgs VEV remains approxi¬ 
mately constant at short times, until reheating is suf¬ 
ficient to destabilize the second minimum, whereas in 
IC-2 the field relaxes as soon as the Hubble parameter 
becomes sufficiently small. Subsequent oscillations have 
a larger amplitude in the IC-1 scenario; this is due to 
the difference in Aj values, which result in less Hubble 
friction, and that the additional term (3) contributes to 
the velocity of the VEV. 

A notable feature in 5 is that shortly after maximum 
reheating, the Higgs condensate begins to oscillate more 
rapidly. This is due to the non-perturbative decay of the 
Higgs condensate, as illustrated in Fig. 3 above. (In the 
IC-2 scenario, such features are not relevant due to the 
rapid decay of the amplitude of oscillation.) 

We see that both the thermal decay of the conden¬ 
sate and the non-perturbative decay of the condensate 
have little effect on first approach of the VEV to zero; 
in the scenario we outline below, the lepton asymmetry 
is generated primarily during this swing, and therefore, 
these processes have little effect on the total asymmetry 
generated. 
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V. EFFECTIVE CHEMICAL POTENTIAL 

In the Introduction, we observed that the Higgs poten¬ 
tial may be sensitive to the effects of higher dimensional 
operators, which are normally suppressed by powers of a 
high scale. In section II, we have seen how such opera¬ 
tors can be used to make a quasistable minimum in the 
Higgs potential or to suppress the growth of the Higgs 
VEV until the end stages of inflation. Now, we consider 
an operator, involving only Standard Model fields, which 
generates an effective external chemical potential for lep¬ 
tons (and also baryons). This operator is 

= (45) 

where is the fermion current of all fermions which 
carry SUl(2) xUy(1) charge. We observe that the zeroth 
component of (45) is the B + L charge density. 

We now consider how an operator of this form can be 
generated. Within the Standard Model itself, one can 
use quark loops and the CP-violating phase of the CKM 
matrix [39, 40] to generate an effective operator of the 
form 

where W and A are the SUl(2) and Uy(1) gauge fields 
respectively. This term is small due to the small Yukawa 
couplings and small CP-violating phase. 

However, a term of the same form can be generated by 
replacing some or all of the quarks with heavier fermions, 
which may have larger Yukawa couplings and/or CP- 
violating phases. The scale in the denominator may be T, 
due to thermal loops, or the mass scale of these fermions, 
Mn [39-42]. In the latter case, it is important that the 
fermions not acquire masses through the Higgs mecha¬ 
nism; otherwise, the Higgs VEV dependence in this term 
cancels out. Such fermions may acquire soft masses sim¬ 
ilarly to higgsinos and gauginos in supersymmetric mod¬ 
els. 

This operator could be generated in a UV-complete 
model. As a concrete example, we mention the fully 
renormalizable Lagrangian 

Chd = 

+ + to'02 V'2 + h.c., (47) 

where i/ii are a set of SU(2) doublets, while ^^2 is a singlet 
under both SU(2) and U(l). Despite the explicit mass 
terms, this Lagrangian is invariant under SU(2) rotations 
provided that both right and left components of the ipu 
doublets couple vectorially to gauge bosons. The phases 
of the tpii doublets may be fixed by eliminating the phases 
in the mass matrix. Provided that * > 3, there are more 
phases 6 i than can be eliminated by rotating the Higgs 
field (j) and the singlet tp 2 - Fermionic loops such as those 
in [39, 40], which involve sums of the Yukawa couplings 


due to insertions of the Higgs VEV (</) generate an 
effective operator of the form (46) . In this case, the scale 
in the Oq operator is A„ ~ M ~ m. 

Once an effective operator of the form (46) is gener¬ 
ated, it may be transformed into (45) through the elec- 
troweak anomaly equation [11]. However, this is only jus¬ 
tified if the electroweak sphalerons are in thermal equi¬ 
librium [43, 44]. Otherwise, the operator (46) involves 
the Chern-Simons number density, which is not changed 
by Higgs relaxation unless the phase of the Higgs VEV 
evolves. At least for slowly evolving Higgs VEVs, the 
sphaleron transition rate per unit volume at finite tem¬ 
perature is 

r,p = exp(-?;/2r), (48) 

where the exponential factor accounts for the suppression 
due to being in the broken phase. As both the Higgs 
VEV and the temperature are quickly evolving in the 
scenario considered here, it may be difficult to arrange 
for the electroweak sphalerons to be in thermal equilib¬ 
rium. However, additional gauge groups which couple to 
fermions can contribute to the anomaly and generate the 
requisite term, as discussed in Appendix A in [12]. For 
our purposes, we simply consider a scenario with opera¬ 
tor (45) without specifying the mechanism by which it is 
generated. 

Returning to equation (45), we observe that integrat¬ 
ing by parts and dropping an unimportant boundary 
term gives 



In the case where A„ = a constant (for example, 
the mass scale of a fermionic loop, as outlined around 
Eq. (47)), this becomes 

= (50) 

If thermal loops generate this term instead, then this 

becomes 

= (51) 

provided that the temperature is slowly varying on the 
time scales of the Higgs oscillation. In the IC-1 scenario 
specifically, the Higgs VEV remains trapped until there 
is sufficient reheating, which generally ensures that the 
temperature will be slowly varying during the evolution 
of the Higgs condensate. Since the Higgs VEV varies 
only in time, these equations become 

06.A„ = — ^(9o/>^)jb+z,. (52) 

For each fermionic species, its contribution to this term 

can be combined with its kinetic energy term, 
which is equivalent to the replacement 

ido ido - {doct?)/ 


(53) 
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FIG. 7. Some diagrams that contribute to lepton number 
violation via exchange of a heavy Majorana neutrino. 


This effective raises the energy of antiparticles, E —> 
E + (docj)'^)/A'^, while lowering it for particles, E —>■ 
E — (i9o0^)/A^. This can be interpretted as an exter¬ 
nal chemical potential; further remarks along these lines 
are discussed in Appendix A. In the presence of a lepton- 
number-violating interaction, the system will relax to its 
equilibrium state, in which the number of particles ex¬ 
ceeds the number of antiparticles. For future reference, 
we define the effective external chemical potential. 


Eo = 



(54) 


As an effective chemical potential, this operator spon¬ 
taneously breaks not only CP, but in fact, CPT [45]. 
This operator has been used previously in spontaneous 
baryogenesis models utilizing gauge [11, 35, 46-48] or 
gravitational [49] interactions. 


VI. LEPTON NUMBER VIOLATING 
PROCESSES 


The universe can relax to its equilibrium state with 
nonzero lepton and baryon number only if there exists 
some lepton-number or baryon-number violating process. 
In order to induce such processes, we consider a minimal 
extension of the Standard Model with the usual seesaw 
mass matrix in the neutrino sector [50] . In theories with a 
nonzero Majorana mass, the effective lepton number L is 
the sum of the lepton numbers of the charged leptons and 
the helicities of the light neutrinos. This is conserved in 
the limits —>• oo and 0, but it is not conserved 

for a finite Mr- Insertions of the Majorana mass induce 
lepton-number-violating processes such as those shown 
in Fig. 7. 

We further require that the Majorana mass be sig¬ 
nificantly greater than both the maximum reheat temper¬ 


ature and the initial mass of the Higgs bosons within the 
condensate (meff((/)o))j which suppresses the production 
of right-handed neutrinos both from thermal production 
and the decay of the Higgs condensate. Consequently, the 
contribution from the typical leptogenesis scenario [51] is 
strongly suppressed. 

The lepton-number-violating diagrams shown in Fig. 7 
necessarily involve the exchange of the heavy right- 
handed neutrino in order to violate lepton number, and 
therefore are comparatively suppressed, leading to a nat¬ 
urally small value for the asymmetry. 

In order to calculate the thermally averaged cross sec¬ 
tion {(jv) for these processes in the early universe, we 
need to know the number densities of neutrinos and Higgs 
bosons. These can be produced directly through the de¬ 
cay of the inflaton, or in the thermal plasma through 
weak interactions, involving weak bosons with masses 
mw oc Generically these weak interactions may 

be in or out of equilibrium the plasma created by infla¬ 
ton decay; however, when ^ 0, these interactions 
will be in equilibrium and equilibrate the distributions 
of charged and neutral leptons. To be concrete, we will 
use a thermal number density of each of these species. 
The calculation of the cross section and reaction rate are 
given in Appendix B. 

We note that /AIp> is set by the mass scale of the 
left-handed neutrinos, such that 


2Mr 


= 0.1 eV. 


(55) 


The cross section found in Appendix B is to a good ap¬ 
proximation a function of y^ jM r only. There is a reso¬ 
nance in the s-channel contribution to the cross section; 
however, we found numerically that this resonance does 
not change the result appreciably. This is not unexpected 
as the energy scale, which is set by the temperature, 
remains significantly below the right-handed Majorana 
mass scale at all times. 

As we will discuss below, sphaleron processes later con¬ 
vert this lepton charge asymmetry into a baryon asym¬ 
metry, as in the typical leptogenesis scenario [51]. 


VII. BOLTZMANN TRANSPORT EQUATION 

The reactions discussed in Section VI are generally not 
sufficient to establish equilibrium, due to the suppression 
from the large Majorana mass. (Recall that we have as¬ 
sumed ^0 ^ and Tmax ^ in order to suppress 
the typical leptogenesis mechanism.) The relaxation of 
the system towards equilibrium can be described by a 
system of Boltzmann equations, based on detailed bal¬ 
ance. The rate of change in the neutrino number density 
is [52] 
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(56) 


where 7®^(7l —)■ i?) is the equilibrium spacetime rate for the process A ^ B. We will assume that interactions are 
sufficiently fast that the Higgs bosons have their equilibrium density, and in equilibrium, the rate for the process 
A ^ B is equal to the rate of i? —H. Therefore this simplifies to 


+ 3i7n., = - ^ 7'=^(z.ih0 77 P,h°) + ( - l) 77 h°h°) 
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(57) 


while for antineutrinos we find the similar equation 


+ 377np^ = - ^ f ^ Y'‘{'^Lh° 77 Vth^) + - 1^ ^‘"‘‘{vLVt 77 /l°/l°) 
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(58) 


Since we are interested in the order of magnitude of the final asymmetry, we simplify to the case in which there is 
only a single neutrino species. Subtracting Eq. (57) from Eq. (58) gives a Boltzmann-type equation for the difference 
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(59) 


The rates 'y^'^{A 77 B) refer to the process A <->■ B 
in equilibrium, but in the presence of the Oq opera¬ 
tor, which alters the energy of particles and antiparticles. 
Consequently, these reaction rates are not generally equal 
to the rates one would find in the absence of the Oq oper¬ 
ator; however, the difference appears at a higher order in 
Eq/T [43] and so we will neglect it. This has the conse¬ 
quence that the rates for 77 i^l^l and hPh^ 77 
are equal. We will use the subscript 0 to denote reaction 
rates calculated without the Oq operator. 

We next substitute n®® = and n?® = 

g-^o/T^eq^ where tt-q® = jn'^ is the equilibrium number 
of left-handed neutrinos (or antineutrinos), when Eq = 0. 
Expanding the resulting equation to lowest order in Eq/ 
T gives 


where we have introduced the notation -frip^, 

and we have dropped terms quadratic in the asymme¬ 
try (e.g., n\). Approximating ^ the equation 

becomes 

flL + SHul = (^riL - ['jQ'^{DLh° 77 l^Lh°) 

+Jo'^ivLi2L ^ h°h°)] . (61) 

The reaction rates are calculated in Appendix B. Erom 
this equation, we observe that the equilibrium asymme¬ 
try is 


2EQ^eq _ 2 T 2 
T ""o 7r2 A2 ■ 


(62) 


flL+SHriL = - ^ 77 ^Lh°) 

- (n/Z^riL - ^ h°h°), (60) 


During subsequent oscillations of the Higgs VEV, the 
chemical potential changes sign. However, due to the 
large suppression in the cross section, significant washout 
can be avoided if the Higgs oscillation amplitude de¬ 
creases rapidly. This is in contrast to Ref. [47], in which 
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expands the analysis of [1], which only considered two 
such scenarios. In all scenarios, we use the improved 
Boltzmann equation (61), with the cross sections calcu¬ 
lated in Appendix B, and the improved calculation of the 
Higgs condensate equation of motion. We show the time- 
evolution of the lepton asymmetry in all four scenarios; 
subsequently, we present an analysis of the parameter 
space in which a sufficiently large late-time lepton asym¬ 
metry can be generated. 

An analytic approximation for the asymmetry calcu¬ 
lated here numerically can be found in [1], which we 
summarize here. The Boltzmann equation (61) can be 
analyzed in two regimes: during the relaxation of the 
Higgs vacuum expectation value, during which is 
significant and Eo{t) ^ 0, and the subsequent cooling of 
the universe, during which Eq = 0. The reactions shown 
in Fig. 7 are typically out of thermal equilibrium by the 
end of Higgs relaxation, due to the exchange of the heavy 
right-handed neutrino, which suppresses washout. One 
may approximate the potential as H ^ with an effec¬ 
tive running coupling A as in [17]. The final asymmetry 
1 = »^L/(27r^5*T^/45) is approximately 


45 -s/A^qA/ 2 t^ 2 ,, ■ fi TiS 1 

V = ^ J^2rp2 X 


X exp 


24 -^ 3^15 


(JrMr\Tr 


(63) 


FIG. 8. The comoving density of equilibrium lepton asym¬ 
metry for IC-1, with the parameters A/ = 10^® GeV and 
F/ = 10® GeV. Purple (Blue) line corresponds to the result 
with (without) the thermalization through the top quark. The 
top diagram corresponds to times before maximum reheating, 
whereas the bottom diagram corresponds to times after max¬ 
imum reheating. 


where frix and Tj-ix as the time and temperature at the 
end of Higgs relaxation, and aR « 10“^^ GeV“^ approx¬ 
imates the the cross section given by equation (B5). This 
estimate includes the dilution due to entropy production 
during the ongoing reheating process. 


washout was avoided by using coherent oscillations of the 
inflaton field to modify the sphaleron transition rate. 

We note that this is modified by the decay of the Higgs 
condensate; however, as discussed above, the Higgs con¬ 
densate does not typically thermalize until after reheat¬ 
ing. Fig. 8 demonstrates the effect of thermalization on 
the equilibrium density. However, since the lepton asym¬ 
metry will be generated primarily during the first oscil¬ 
lation of the Higgs VEV, the effect of the thermalization 
of the Higgs condensate is negligible. 


VIII. RESULTING ASYMMETRY 

In this section, we consider the lepton asymmetry pro¬ 
duced by these Higgs-neutrino interactions during the 
relaxation of the Higgs VEV, as outlined above. We 
present four numeric examples, covering both IC-1 and 
IC-2, along with the scale of the operator A„ set to the 
temperature T (motivated by thermal loops) and a con¬ 
stant Mn (motivated by loops of heavy fermions). This 


A. Four Numerical Examples 

In this subsection, we present the lepton asymmetry as 
a function of time for the four scenarios mentioned above. 
First, we consider two scenarios for IC-1 in Fig. 9, one 
with A„ = T (blue, solid) and one with A„ = = 10^^ 

(red, dashed) for the relevant scales in the Oq opera¬ 
tor. Both scenarios have a maximum temperature of 
6 X 10^^ CeV, since they share the inflationary param¬ 
eters A/ = IQis CeV and F/ = 10® CeV. As in Fig. 5, 
the initial Higgs VEV is 10^^ CeV in both cases, which is 
set by the location of the second minimum in the Higgs 
potential. Although the asymmetry rj oscillates during 
the first few oscillations of the Higgs VEV, it relatively 
quickly settles into a steady state, and approaches a con¬ 
stant value around the beginning of the radiation dom¬ 
inated era. Note that the Higgs field begins to oscillate 
before the time of maximum reheating. 

As mentioned above, the cross section depends primar¬ 
ily on y'^/M r which is fixed by the light neutrino masses. 
As mentioned above, we require T <C Mr in order to sup¬ 
press the thermal production of right-handed neutrinos; 
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FIG. 9. Plot of the resulting asymmetry for IC-1, for A„ = 
T (blue, solid) and An = Mn = 10^"* GeV (red, dashed). 
Both scenarios have A/ = 10^® GeV, and Tj — 10® GeV. The 
vertical lines designate the first time the Higgs VEV crosses 
zero, time of maximum reheating, and the beginning of the 
radiation dominated era, from left to right, t = 0 corresponds 
to the beginning of inflaton oscillations. 




FIG. 10. Plot of the resulting asymmetry for IC-2, for A„ = T 
(blue, solid) and An = M„ = 5 x 10^® GeV (red, dashed). 
Both scenarios have A/ = 10^® GeV and Pj = 10® GeV. 
From left to right, the dotted lines correspond to the time 
of maximum reheating, the first time the Higgs VEV crosses 
zero, and the beginning of the radiation dominated era. 


FIG. 11. The resulting asymmetry (log|? 7 |) at the end of 
reheating for IG-1, for An = M„, with Aj = 10^® GeV. 


to determine the Higgs VEV at the end of inflation; this 
results in (j)Q = 10^^ GeV for the Higgs VEV as the start 
of Higgs relaxation. (We emphasize that this choice, with 
Mn < 4>o and Mn < T, raises questions regarding the use 
of effective field theory, which we address below.) 

In order to suppress the thermal production of right- 
handed neutrinos, we have taken Mr = lOTmax = 
10^® GeV; in order to produce left-handed neutrino 
masses on the scale of 0.1 eV, the neutrino Yukawa cou¬ 
pling must be 1.9. (This gives y®/47r « 0.3.) 

The final asymmetries here are of order 10“^^ (for 
A„ = T) and IQ-^^ (foj. A„ = M„ = 5 x 10^® GeV). 
As Mn is generally smaller than the temperature, it is 
not surprising that this results in a larger asymmetry. 
These values are insufficient to account for the observed 
matter-antimatter asymmetry; this motivates a search of 
the available parameter space. 


we found that it was sufficient to set Mr = 9 x 10^® GeV, 
which results in ~ 1.7 using equation (55). (This gives 
jAm ^ 0.2, within the perturbative regime.) 

The late time asymptotic asymmetry is rj ^ 10“^ for 
An = T and ry ~ 10“® for A„ = = 10^'^ GeV; this is 

expected as the temperature is lower than M„. We dis¬ 
cuss the variation of the final asymmetry over parameter 
space below. 

First, however, we present similar results for the IG-2 
scenario, again for the two cases A„ = T (blue, solid) 
and An = Mn = 5 x 10^® GeV (red, dashed). Both plots 
have the inflationary parameters Aj = 10^® GeV and 
T/ = 10® GeV, which results in a maximum temperature 
of 10^"^ GeV during reheating. We again take Mast = 8 


B. Parameter Space 

In two of the four scenarios above, the resulting lepton 
asymmetry is 0(10“®) or larger, which is sufficient to 
explain the observed baryon asymmetry. However, it is 
interesting to explore the resulting asymmetry as a func¬ 
tion of parameter space; results are shown in Figures 11, 
12, 13, and 14. 

As above, we handle the initial conditions with the op¬ 
erator and scale given in HA for the IC-1 plots, and as 
discussed in HB with Mast = 8 for the IC-2 plots. We 
emphasize again that the resulting asymmetry is sen¬ 
sitive to u'^/Mr, which is set by the left-handed neu- 
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FIG. 12. The resulting asymmetry (log|? 7 |) at the end of 
reheating for IC-2, for A„ = with A/ = 5 x 10^® GeV, 
which gives 0o = 2.7 x 10^“^ GeV. 


trino mass scale, and not on the specific value of Mr. 
However, to suppress thermal production of right-handed 
neutrinos, we have chosen Mr = lOTmax (for IC-1) and 
Mr = 20Tniax (for IC-2). We have then set the neutrino 
Yukawa coupling y by the scale of the left-handed neu¬ 
trino masses (equation (55)). We have noted in gray the 
regions in which the perturbativity condition j/^/47r < 1 
fails. 

For the IC-1 plots, the post-infiationary Higgs VEV 
00 is determined entirely by the operator which lifts the 
second minimum to generate the quasistable vacuum; for 
the operator and scale in HA, the Higgs VEV relaxes 
from 00 = 10^® GeV. For IC-2, 0o is determined by the 
Hubble parameter during inflation, which is in turn fixed 
by the energy density in the inflation field (see equation 
( 6 )); 

First, we remark on some general features. The asym¬ 
metries generated in the IC-2 scenario are smaller than 
those generated in the IC-1 scenario. This is because in 
IC-1, the Higgs VEV does not evolve until the tempera¬ 
ture is sufficiently large to destabilize the false vacuum; 
consequently, the initial evolution of the VEV to zero oc¬ 
curs at higher temperatures. (Compare the vertical lines 
significantly the first Higgs VEV crossing and maximum 
reheating in Figures 9 and 10.) As a result of the higher 
temperature, the system is driven towards equilibrium at 
a faster rate (through the Boltzmann equation (61)); fur¬ 
thermore, in the A„ = T scenario, the larger temperature 
also means that the equilibrium charge density is larger. 

Figures 11 and 12 show the lepton asymmetry ry as a 
function of parameter space, in the case in which the scale 


of the Oq operator is a constant M„. To reach compara¬ 
ble asymmetries in the IC-2 scenario, we must decrease 
the scale Af„ significantly, such that throughout this plot, 
Mn < 00 and M„ < Tmax- In the IC-1 plot, these con¬ 
ditions fail below the red dashed line and blue solid line 
respectively. In these regions, the use of effective field 
theory in generating the operator (45) is questionable. 
An ultraviolet completion of the model is necessary to ob¬ 
tain a reliable description of the dynamics in the regime 
where the temperature exceeds the scale M„. We leave 
such a completion, which would also elucidate the nature 
of the new physics leading to the appearance of the Oq 
operator, for a future work. 

We focus on the region of Fig. 11 for which the asym¬ 
metry r] is larger than 10“^° and Mn > O.lTmax- We 
see that this favors smaller values of F/. However, 
for a given A/, there is a minimum F/, for which the 
maximum temperature is insufficient to destabilize the 
second vacuum. For the parameters considered here 
(A/ = 10^® GeV and the lift operator given in HA), 
this occurs for F/ = 6.3 x 10®. 

Next, we consider the case in which the scale of the Og 
operator is set by the temperature, in Figures 13 and 14. 
This parameter space has one fewer parameters, and so 
we allow A/ to also vary, which changes the Hubble pa¬ 
rameter during inflation. For IG-2, increasing Hj results 
in a larger value of 0o, as described by (6), which in¬ 
creases the resulting asymmetry. This also increases the 
temperature scale, resulting in a larger asymmetry, as is 
evident in both figures. (We also note that for IC-1, we 
must take care that quantum fluctuations during infla¬ 
tion do not destabilize the second vacuum; this is shown 
in orange in Fig. 13.) 

As mentioned above, if the reheat temperature is suf¬ 
ficiently small, thermal corrections are unable to desta¬ 
bilize the second vacuum, and therefore this is no relax¬ 
ation of the Higgs VEV. This region is denoted in white 
in Fig. 13. Furthermore, in the region in which Mr < 0o, 
right-handed neutrinos can be copiously produced by the 
decay of the Higgs condensate, which is not desirable (as 
concerns the lepton asymmetry production scenario pre¬ 
sented here); this region is denoted in yellow. Further¬ 
more, if A/ is too small, there is insufficient inflation to 
account for the observed flatness and uniformity of the 
universe; this region is shown in blue on both figures. 

In IC-2, there is a further concern that the Higgs VEV 
can probe the second, deeper minimum at large VEVs. 
This may not be a phenomenological problem [53], but 
would require a refinement of the analysis presented here. 
(Alternatively, Mast could be decreased, such that 0o 
remains below the instability scale.) This region is shown 
in purple in Fig. 14. 

We see that for IC-1, it is possible to find parameter 
space in which a sufficiently large asymmetry is gener¬ 
ated, but this is not possible for IC-2. For IC-1, smaller 
F/ values are favored (and consequently, slower reheat¬ 
ing), as for constant A„. 
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of freedom go out of thermal equilibrium; the resulting 
entropy production reduces the asymmetry by about two 
orders of magnitude. 

This process has produced a net density of {B — L) 
charge, which is unchanged once these processes are neg¬ 
ligible. However, the {B + L) C/(l) symmetry is anoma¬ 
lous, and electroweak sphalerons will redistribute the ex¬ 
cess between leptons and baryons as in standard lepto- 
genesis [51], at a rate per unit volume 

Tsp ~ (awT)^ exp [-gw(l)it)/T]. (64) 

At small vacuum expectation values, the B and L 
densities approach their equilibrium values, ng = 
{28/79)nB-L- This produces a baryon asymmetry of 
about the same order of magnitude as the lepton asym¬ 
metry found above. Consequently, the regions of param¬ 
eter space that generate ij ^ 10“® in the analysis above 
give a final baryon asymmetry matching the observed 
value of 0(10-1°). 


FIG. 13. The resulting asymmetry (log|? 7 |) at the end of 
reheating for IC-1, for A„ = T. 


> 

o 

< 

'5o 

o 


log(r, [GeV]) 

FIG. 14. The resulting asymmetry (log|? 7 |) at the end of 
reheating for IG-2, for A„ = T. 


C. Converting the Lepton Asymmetry Into a 
Baryon Asymmetry 

Thus far, we have analyzed the production of an excess 
of leptons over antileptons; here we discuss how this is 
converted into a baryon asymmetry. First, though as the 
universe continues to cool, the Standard Model degrees 


IX. CONCLUSIONS 

In this paper, we have extended the analysis of Ref. [1], 
which introduced a novel leptogenesis possibility in which 
the lepton asymmetry is as a consequence of an effective 
chemical potential induced by the post-inflationary relax¬ 
ation of the Higgs field. Although right-handed neutrinos 
participate in the lepton-number-violating interactions as 
a mediator, this is different from the typical leptogene¬ 
sis scenario in which the asymmetry is produced via the 
decay of right-handed neutrinos. Even though the heavy 
right-handed neutrino suppresses the cross section which 
produces the asymmetry, we have shown parameters for 
which a sufficiently large the asymmetry is generated. 

We have analyzed the evolution of the Higgs conden¬ 
sate in detail, including both non-perturbative and per¬ 
turbative decay. We have derived the relevant Boltzmann 
equation which governs lepton number, and we have re¬ 
placed the order-of-magnitude estimate with a tree-level 
scattering cross section between Higgs bosons and neu¬ 
trinos in the thermal plasma. Furthermore, we have con¬ 
sidered the evolution of the lepton asymmetry for four 
combinations of producing the large Higgs VEV during 
inflation (IC-1 and IC-2) and the scale of the Oq opera¬ 
tor (a fermion mass scale Mn and the temperature T); we 
then presented an analysis of the asymmetry as a func¬ 
tion of parameter space. We demonstrated regions which 
produces a baryonic asymmetry that meets or exceeds 
observational limits. 
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Appendix A: Interpretting the CJg Operator as an 
External Chemical Potential 

In section V, we remarked that the Og operator in 
equation (46) acts like an external chemical potential. In 
this appendix, we explain why this is so and how this 
leads to a number density asymmetry in chemical equi¬ 
librium. 

This Og operator induces a term proportional to 
{do(f^)lK,3B+L in the Lagrangian. If (j) is treated as 
an external field (which we discuss further below), then 
this produces a term of the form in the 

Hamiltonian, which has the appropriate form —/ieffJs+i- 

A term similar to this, using the phase of the Higgs 
VEV is frequently used in spontaneous baryogenesis sce¬ 
narios, in which the phase of the Higgs VEV is used in¬ 
stead of its magnitude (e.g., [54]), 

o; = (Ai) 

However, in such scenarios, the asymmetry is produced 
via the decay of the Higgs condensate, and therefore, it 
is not appropriate to treat 9 as an external degree of 
freedom. When the Hamiltonian is determined using 


space densities are 

fp = exp(-(£; -Eg- fJ.p)/T) 
fp=exp{-{E + Eg-fj.p)/T) (A3) 

The number densities of particles and antiparticles can 
be found in the normal manner, using 

^P = J exp{-{E -Eg- Hp)/T) 

np = j-^^^^V{-{E + Eg-p.p)/T) (A4) 

If we use the non-relativistic relation E = I2m^ then 
we find 

fip = -Eg -I- Tln(A^np) 

Pp = Eg+ T In(A^np), (A5) 

where A = \/27rmT. In the above relation, the first term 
can be interpretted as an external chemical potential (due 
to the “driving” effect of the Og operator), while the 
Tln(A^np) is the usual chemical potential of an ideal gas. 

If a lepton-number-violating process or baryon- 
number-violating establishes chemical equilibrium be¬ 
tween the species, then the chemical potentials will be 
equal, Pp = Pp- This gives the expected result 


A similar result can be derived using the relativistic re¬ 
lation E = p instead. 




(A2) 


there is no contribution from Og. Although an asym¬ 
metry may be produced in such cases [55-57], it is not 
appropriate to interpret d as a chemical potential. 

In the scenario we consider in this work, the time scale 
for the reactions which maintain the thermal distribution 
of the plasma is smaller than that of the evolution of 
the Higgs VEV. Therefore, for purposes of asymmetry 
generation, it is reasonable to consider the Higgs VEV 
as a background field, in which case it is appropriate to 
consider this as a chemical-potential-like term [57] , as we 
explain below. 

The Og operator shifts idg —>■ idg — {dg(f)^)/K^ in the 
Lagrangian. Consequently, the asymptotically free eigen¬ 
functions are ^ exp(=Fj(E=F (i9o(/'^)/A^)t), which justifies 
our comment that this is equivalent to decreasing the en¬ 
ergy of particles by Eg = {dg(j)^)/K\ and increasing the 
energy of antiparticles by the same amount. 

If we use the ideal gas approximation, then the phase 


Appendix B: Calculation of 
Lepton-Number-Violating Cross Section and 
Reaction Rate 

In this section, we calculate the cross section and re¬ 
action rate for the processes shown in Fig. 7, assuming 
a thermal number density for Higgs bosons and neutri¬ 
nos, as discussed in Section VI. This improves the order 
of magnitude estimates used in [1]. As explained in the 
text, we can use the approximate cross section with the 
energy shift due to the Og operator set equal to zero. In 
this approximation, the reaction rates for processes with 
neutrinos and antineutrinos are equal. 

The top two diagrams of Fig. 7 are the s- and t-channel 
diagrams of the process /i°iz —)• whereas the bottom 

diagram describes the process lyv hPh^. The s-channel 
has a resonance at E ~ Mr] however, the typical energy 
scale T is far beneath this. For completeness, we include 
the resonance, although it will have negligible effect. In 
calculating these cross sections, we follow the conventions 
of [58] for the Feynman rules of Majorana fermions. 

The matrix element for the process is 
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—iM = i 


YL^Y* 


Mm — iTi/2 
s — -f 




Mm — iTi/2 


XLa(pi,Si)y'l{pi,Si)5p, (Bl) 


where s and t are the Mandelstam variables, and is 
the width of the right-handed Majorana neutrino. (For 
a discussion of Breit-Wignar propagators, see [59]). The 
indices 1, 2, 3, and 4 refer to the incoming neutrino, in¬ 
coming Higgs boson, outgoing Higgs boson, and outgoing 
antineutrino, in that order. The index i indicates a sum 
over the heavy right-handed Majorana neutrinos. Let us 
define 


Ai — s — M/ii + r^/4, 

Bi = t-M/,^ + r^ji, 

C, = T,Mm- (B2) 


Then the matrix element squared, summed over both 
the initial and final spin states (as discussed in [52]), is 


y:iA<i“ = y: 2 p. 

Si,S2 i 


X P4 


\Ym?\Yu? 




2{A,B, + Cf) 


A^ + Cf Bf + Cf {A^Bi + Cfy + CfiAY 


X 


B^r\ 

(B3) 


In the center of mass reference frame, the cross section 
is 


-om{s)=—y: 


\Ym?\Yu? 


Ml, 


dt (s -I- t) 


1 


1 


,A^ + Cf Bf + Cf 


2{A,B, + Cf 
{A,B, + C/Y + C/{A, - BiY\ ■ 

(B4) 


Generically, the thermally averaged cross section is related to the CM cross section by 


(cru) = 


_ \ _ /■“ [s - (mi - m 2 f][s - (mi -H mz)^] /rc', 

Q'T' 2 7^ / /'T'\ X/ 2 / l'T'\ I r~ v 

8T X mfK 2 {mi/T) x m^K 2 {m 2 /T) VS 


and so the thermally averaged cross section for h^iy —> h^v is 


{(j{h^v h°9)v) = 

i 

1 


\Ym?\Yu? f^2^ ^ C 


5127r 


dx dy{x^ - y^Kiix) 


(x 2 T 2 -Ml, + rf/4)2+rfM|, 


2[{x^T^ - Ml, + rf/4)(y^r^ + - r|/4) - 


\y^T^ + - rf/4)2 + ((^ 2^2 _ Ml^ + rf/4)(y2T2 + - Vj/A) - VjMl^Y + T^MUx^ + y^YT\ 

fB6) 


where we have introduced the dimensionless variables 
X = \fs/T and y = Y^A/T. Since the temperature 
evolves in time, the cross section also does; however, 
when expanded in powers of T/Mm, the lowest order 
contribution is ~ 1/M^j, as expected. Repeating the 
same steps with the I'ii'L —> cross section, which 

does not have a resonance, gives 


{u{viVL h°h°)v) = y] 

i 


\Ym\^\Yu\^ 

647rMp- 


(B7) 


The reaction rates are related to these cross sections 
by 


Y'^{al3 yd) = {rTY){n/Y) (cr(a/3 -i5)v ), (B8) 


which holds for any 2 —)■ 2 process. Since we take ifo = 0 
in this section, the number densities for Higgs bosons, 
neutrinos, and antineutrinos are all equal to and 

so for both processes. 


7 


0 



(B9) 


As noted in the text, in order to simplify the calcula¬ 
tion, we will consider only the case in which the flavor 
indices ^ and L are equal, and the contribution of a sin¬ 
gle right-handed neutrino dominates. Its decay rate is 
r ~ y’^Mji/lQ-n, from the only decay Nfi h^VL- 
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